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ABSTRACT 

For gamma-ray bursts (GRBs) born in a stellar wind, as the reverse shock crosses the ejecta, usually 
the shocked regions are still precipitated by the prompt MeV 7— ray emission. Because of the tight 
overlapping of the MeV photon flow with the shocked regions, the optical depth for the GeV photons 
produced in the shocks is very large. These high energy photons are absorbed by the MeV photon 
flow and generate relativistic pairs. These pairs re-scatter the soft X-ray photons from the forward 
shock as well as the prompt 7— ray photons and power detectable high energy emission, significant 
part of which is in the sub-GeV energy range. Since the total energy contained in the forward shock 
region and the reverse shock region are comparable, the predicted sub-GeV emission is independent 
on whether the GRB ejecta are magnetized (in which case the reverse shock IC and synchrotron self- 
Compton emission is suppressed). As a result, a sub-GeV flash is a generic signature for the GRB 
wind model, and it should be typically detectable by the future Gamma-Ray Large Area Telescope 
(GLAST). Overlapping also influence neutrino emission. Besides the 10 15 ~ 10 17 eV neutrino emission 
powered by the interaction of the shock accelerated protons with the synchrotron photons in both the 
forward and reverse shock regions, there comes another 10 14 eV neutrino emission component powered 
by protons interacting with the MeV photon flow. This last component has a similar spectrum to the 
one generated in the internal shock phase, but the typical energy is slightly lower. 

Subject headings: gamma rays: bursts - acceleration of particles - elementary particles 



1. INTRODUCTION 

The leading model for long gamma-ray bursts (GRBs) 
invokes a relativistic jet emerging from a collapsar (e.g. 
Woosley 1993; Paczynski 1998). It has been widely ex- 
pected that GRBs associated with supernovae should oc- 
cur in a pre-burst stellar wind enviroment with particle 
number density n oc R~ 2 (e.g. Dai k Lu 1998; Meszaros, 
Rees & Wijers 1998; Chevalier & Li 2000). As the jet 
expands into a dense steller wind, a forward shock (FS) 
and a reverse shock (RS) form (Chevalier & Li 2000; Dai 
& Lu 2001, hereafter DL01), whose synchrotron emission 
peaks in the soft X-ray band and in the ultraviolet to op- 
tical band, respectively. Photomeson interaction in the 
RS region also results in neutrino emission in the energy 
range of 3 x 10 15 eV - 3 x 10 17 eV (DL01). 

However, until very recently, in nearly all the works on 
RS emission, possible overlapping of the initial prompt 
7— rays emission (MeV photon flow) with the RS region 
has been ignored. Beloborodov (2005) pointed out that 
when overlapping is important, the RS emission would 
be suppressed significantly since the electron cooling is 
dominated by inverse Compton (IC) scattering off the 
MeV photon flow. Correspondingly, the upscattering 
would power a strong GeV-TeV photon flash. Fan et 
al. (2005a) suggested that when overlapping is substan- 
tial, similar processes are also relevant to the early FS 
emission. 

In this work, we show that overlapping is a common 
feature for GRBs born in a steller wind (SJ21 7 so that the 
early afterglow wind model needs substantial revision. 
We then systematically study the consequence of this 



overlapping effect in the wind model. In fJ3we show that 
such an early photon-shock interaction inevitably leads 
to the prediction of a sub-GeV photon flash in the wind 
model, which is generally detectable by the Gamma-Ray 
Large Area Telescope (GLAST). In 21 we discuss vari- 
ous high energy neutrino emission processes in the early 
afterglow stage as the result of photon-shock interaction. 
We summarize our conclusions in fj^l with some discus- 
sions. 

2. TIGHT OVERLAPPING OF THE MEV PHOTON FLOW 
WITH THE SHOCKED REGIONS 

The overlapping of the MeV photon flow with the fire- 
ball ejecta is important if the RS crossing radius satisfies 

R x « 2r 2 A, (1) 

where To is the initial Lorentz factor of the ejecta (at 
the internal shock phase), A = cTgo/(l + z), Tgo ~ 20s is 
the observed duration of the long GRB, and z ~ 1 is the 
redshift. The shock crossing radius R x can be written in 
a general form 



R x = max(i? 7 ,r 2 A), 



(2) 



where i? 7 is the radius where the mass of the medium 
collected by the fireball is equal to 1/Tq of the fire- 
ball mass (which corresponds to T ~ To/2), and T 
is the Lorentz factor of the shocked ejecta. In equa- 
tion the first term dominates if the RS is non- 
relativistic, while the second term dominates if the RS 
is relativistic. For the ISM case with a typical number 
density nrgM ~ lcm~ 3 , the RS is non-relativistic and 

B.B X 1U £/ 53 6 n ISM l 2 5 , 



one has R x (ISM) = R n 
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which is usually larger than (or at least comparable 
to) 2I^A = 5.4 x 10 16 cm r 2 , 2 5 Aii. 5 . Throughout 
this paper, we adopt the convention Q x — <5/10 x to 
express physical parameters in cgs units. So, in the 
ISM case, the initial 7— rays can not change the RS 
emission significantly. In the stellar wind case, on the 
other hand, one has n = 3.0 x 10 35 A*i? -2 cm -3 , where 
A* = (M/W- 5 M Q yr^Vw/lO^m s' 1 )" 1 , M is the 
mass loss rate of the progenitor, and v w is the wind ve- 
locity (Chevalier & Li 2000). With such a dense medium, 
the RS is relativistic (Chevalier & Li 2000; DL01): 

F « C 1/4 ro /2 /\/2, where C « 2600^ 53 . 6 /(A* Aii. 5 r^ 2 . 5 ) 
is the ratio of proper mass density of the unshocked 
ejecta and the medium density, E ~ 4 x 10 53 ergs is 
the isotropic energy of the ejecta, which is calculated 
by adopting T90 = 20s, z = 1, and by assuming that 
the observed 7— ray luminosity is L 1 ~ 10 52 ergs s _1 
and that the radiation efficiency of the internal shocks 
is r\ ~ 0.2. The total wind luminosity is therefore 
Ltot — L^/r] = 5xl0 52 ergs s _1 , and the afterglow kinetic 
energy luminosity is L ~ £tot(l — 17) = 4 x 10 52 ergs s _1 
in the beginning of the afterglow phase. The resulting 
bulk Lorentz factor in the shocked region is therefore 

r^87Ll{ 4 6 A: 1/4 . (3) 
The RS crosses the ejecta at the radius 

R x (wind) ~ 2.7 x 10 15 cm£^ 2 6 A* 1/2 A U . 5 . (4) 

At this radius, the rear of the MeV photon flow only leads 
the rear of the ejecta by a distance of R x (wind)/2FQ -C 
A, so that during most of the time when RS crosses 
the ejecta, the shocked regions are still precipitated by 
the MeV photons, so that the electrons and pairs 
in the shocks lose energy mainly by scattering with the 
MeV photons. As a result, the conventional synchrotron 
ultraviolet-optical emission of the electrons in the RS re- 
gion would be greatly suppressed. Let us define a pa- 
rameter Y = Uj/Ub, where U 1 ~ L 7 /47ri? 2 F 2 c is the 
MeV photon energy density in the rest frame of the 
shocked region, and Ub = B 2 /8n ~ (F /r) 2 eBn4m p c 2 
is the comoving magnetic energy density in the same re- 
gion, where m = (1 — rj)L / '4TrR 2 rlm p c 3 is the comoving 
density of the un-shocked ejecta. As usual, we intro- 
duce magnetic field and electron equipartition parame- 
ters £b ~ 0.001 — 0.1 and e ~ 0.1 in the shocked region. 
The wide range of cb is adopted with the consideration 
of broadband afterglow fits (e.g. Panaitescu & Kumar 
2001) and the possibility that the ejecta may be magne- 
tized (e.g. Fan et al. 2002; Zhang et al. 2003). In this 
work, we normalize our expression by taking cb ~ 0.01. 
After some simple algebra, we have 

Yttr)/[(l-r))e B }. (5) 

For e B ~ 0.001 - 0.1 and -q « 0.2, we have Y ~ 
2.5 — 250 ^> 1. The cooling frequency v c (see Eq. (15) 
of DL01) would then be lowered by a factor 1/(1 + Y) 2 . 
For ee,-2 ~ 0.1 — 1, the flux of the soft ultraviolet-optical 
photons is significantly lower than the one calculated in 
DL01, so that the energy loss fraction for each proton 
through pion production (i.e. f n , defined as the ratio 
between the comoving expansion time t' d and the comov- 
ing proton cooling time t' % due to photomeson cooling) 
should be also lowered by roughly a factor of 1/(1 + Y) 



(see equations [21] and [22] in DL01). Consequently, the 
resulting detectable neutrino flux should be somewhat 
lower. 

It is worth pointing out that for some other processes 
accompanying prompt 7— ray emission, such as the pre- 
diction of the ultraviolet flash powered by neutron-rich 
internal shocks (Fan & Wei 2004), the MeV photon flow 
also plays an imporatnt role on the cooling of the ac- 
celerated electrons. By taking into account the over- 
lapping effect, Fan et al (2005b) have shown that with 
reasonable parameters, the near IR, optical and UV emis- 
sion flux from the neutron-rich internal shocks can be as 
strong as lOmJy. Comparing with the standard inter- 
nal shock models, it also suffers weaker synchrotron-self- 
absorption and seems to be able to account for prompt 
long-wavelength observation of GRB 041219a (e.g., Blake 
et al. 2005; Vestrand et al. 2005). 

3. SUB-GEV PHOTON FLASH 

Bcloborodov (2005) suggested that electrons acceler- 
ated by the RS would scatter the MeV photon flow and 
power a detectable GeV-TeV flash. It is suggested that 
this mechanism may potentially acount for the distinct 
high energy spectral component observed in GRB 941017 
(Gonzalez et al. 2003). 

One uncertainty is that, if in the RS region e B > 0.1, 
as might have been the case for GRB990123 (Fan et al. 
2002; Zhang et al. 2003), the Y parameter (equation 
H3) is not large, and the IC effect may not be strong. 
(Hereafter the superscripts r and / represent the RS and 
FS, respectively. If the parameters in both regions are 
the same, no superscript is marked.) Below we will show 
that by taking into account the photon-shock interaction 
in the FS region, the IC process is generic regardless of 
the magnetization of the ejecta. Also, the typical energy 
of the emerged spectrum peaks in the sub-GeV range 
(see below). 

At i? x (wind), the energy contained in the RS region 
can be written as 

E\ ot a j 34 TM ci c 2 « r M ojC 2 /2 « E { tot , (6) 

where we have taken the Lorentz factor of the shocked 
ejecta relative to the unshocked ejecta as 734 ss (ro/r + 
r/r )/2 ~ r /2r > 1 (which is valid for the wind RS 
model). Here M e ; ~ E/{Tqc 2 ) is the rest mass of the 
ejecta. Equation |5J| suggests that the energy contained 
in the RS region is roughly equal to that in the FS re- 
gion. As a result, if a significant part of the FS energy 
is emitted in the sub-GeV range (we'll show below that 
it is the case), whether the sub-GeV emission from the 
RS region is suppressed or not is observationally unim- 
portant. In view of this fact, below we adopt the same 
physical parameters (e c ~ 0.1 and ee ~ 0.01) for both 
the RS and the FS regions. 

The total thermal energies in the RS and FS regions 
are E\ h ~ ( 734 - lKot/734 and £& ~ (r - l)£ t f ot /r, 
respectively. The ratio of them is E\ h / E^ ~ (734 — 
l)/734- For typical parameters, 734 ~ 2, we have £^ h : 
E^ h ~ 1 : 2. This result suggests that for both the high 
energy photon emission and the neutrino emission (as 
discussed in §4), the component from the FS region is 
the dominant one. 

As usual, we assume that the shocked electrons dis- 
tribute as dNe/d'je oc 7~ p (p ~ 2.2) for 7, n < 7 C < 7m, 
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where 7m ~ lO^/B 1 ' 2 is the maximum Lorentz factor of 
the electrons accelerated by shocks, B is the magnetic 
field strength of the shock regions, which can be calcu- 
lated by 



B « 1. 



x 10 3 G e]( 2 _ 2 



A 



-1/4 T 1/4 „_i 



(7) 



For the FS and RS, j m can be estimated by 
7^ = (r - l)e e (m p /m e )[(p - 2)/(p - 1)] 

-1/4 r 1/4 



:2600e e 7 L 



52.6' 



7m = (734 - l)e e (m p /m e )[(p - 2)/(p - 1)] 



•23e n 



4 l/4 7 -l/4 r 



(8) 



^52.6 1 0,2.5- 

Performing IC correction on equation (14) of DL01, we 
also get the cooling Lorentz factor j c < 1, which means 
that the shocked electrons all cool so rapidly that es- 
sentially all of them have 7 C ~ 1. In other words, all 
electrons are in the fast cooling phase. 

The electrons in the shocked regions scatter with the 
MeV photons and produce high energy photons. Given 
the observed prompt GRB peak energy e 7 obs = 300keV 
(throughout the paper, the subscript "obs" denotes the 
parameters measured in the observer frame), the elec- 
tron energy above which the IC process is suppressed by 
the Klcin-Nishina limit is given by (see also Beloborodov 
2005) 

rm e c 2 

In the RS region, one has 7* > 7 m . The peak energy of 
the upscattered photons can be then estimated as 



r,IC 



r 2 b 



1, obs ~ 7m £ 7,obs 



-0.16GcV 



l/2 r 2 



A 1 V 



LI 



0,2.5-^52.6 



1/2/ fc 7,obs 



.(10) 



v 300keV 

In the FS region, on the other hand, one has 7 m 3> 7^, 
so that the IC cooling by the MeV photons is suppressed 
significantly. The fresh energetic electrons lose energy 
mainly by synchrotron self-Compton (SSC) effects, the 
corresponding Compton parameter can be estimated by 
Y$sc — \AcAb ^ 1. The typical energy of SSC photons 
reads 



f.ssc 



^( 7m ) 4 erB/[27rm c(l + z)] 

: 86GeV (1 + «) _1 4-i c B[-a^ 3/a 

r 3/2 R _i 
lj 52.6 rl 15 ' 



(11) 



where h is the Planck's constant. 

Apparently, equations (fTT)ft and fTT|l are quite different. 
However, as shown below, photons more energetic than 
~ 1 GeV will be absorbed by the MeV photons near the 
shock crossing radius, and most energy will reprocessed 
and end up in the sub-GeV range. This is valid for both 
the RS and FS. The pair production optical depth for 
photons with energy ~ 1 GeV (which will be absorbed 
by the soft photons with energy e a ,obs — 2(rm c c 2 ) 2 / (1 + 
z) 2 GeV ~ lMeV) is roughly (e.g. Svensson 1987) 

ll(T T N >e . 



f (lGeV): 



720irR 2 

< 20^,53^ 



)L 



1/2 aV 2 r- 2 



52.6 



15 



(12) 



where iV> Ca oba is the total number of photons satisy- 
ing e bs > £a,obs, which can be estimated by A >(Saob3 w 
0.2J5 7 /[(1 + z)e^ ohs ] for e a:0bs > e 7 , obs , where the coef- 
ficient 0.2 is a rough spectrum correction factor on the 
photon number calculation derived from the integration 
of the GRB spectrum. In this paper, we take a typical 
broken power law GRB spectrum for MeV photons, i.e., 
n(e 7 ) oc e^ 1 for e 7 /50 < e 7 < e 7 , and n(e 7 ) oc e~ 2 for 
e 7 < e 7 < 50e 7 . 

With typical parameters, we have r 77 (lGeV) ~ 2 at 
R x (wind). This means that photons more energetic than 
lGeV are absored by the MeV photons and produce 
pairs. The random Lorentz factor of these secondary 
pairs can be estimated as j± ~ (1 + z)e 7j0 bs/2rTO c c 2 = 
11(1 + z)(e 7 , obs /lGeV). For e 7 , obs = e f ^ ~ 32(1 + 

z) _1 GeV, we have 7^i ssc ~ 350. These pairs will scat- 
ter both the MeV photons and the X-ray photons pro- 

f SSC 

duced in the forward shock region. Since 7^ is sig- 
nificantly larger than 7*, so that the Klein-Nishina cor- 
rection is still important for the pair— MeV photon IC 
scattering. The cross section can be approximated as 
a ps A(x)<jt, where x ~ 7± SSC /7* ~ 7^ ssc /r s» 4, 
and A(x) w ^(hi2x + |) ps 0.2 1 . The energy density 
of the soft X-rays (FS synchrotron emission) is about 
g = (1 - r))e e /[ri(l + Y ssc )} ~ 0.1 times of that of the 
MeV photons, so the fraction of energy of the secondary 
pairs that goes to scattering with the soft X-ray FS emis- 
sion (with a typical frequency defined by eq. |15j ) can be 
estimated as [%/A(4)]/[l + Ag/A(A)] ~ 2/3. The typical 
energy of the upscattered soft X-rays can be estimated 
as (at i? x (wind)) 

f.SSC 

,f,IC f,SSCx2 f n onn„v/^± ^2 



■l7± ) 4, 



obs 



0.29GeV(- 



350 



-y 



1 + Z \-l r 2 ,1/2 
2 ' c,-l fc B,-2 



A* 1 L 52 .q. 



(13) 



The pairs also scatter the initial MeV photons. How- 
ever, the upscattered photons are still too energetic 
to escape. A pair photon cascade likely develops un- 
til the photon energy is low enough to be transpar- 
ent. For example, in the extreme Klein-Nishina limit, 
the typical energy of the upscattered MeV photons 
can be estimated by (1 + z)e obs ~ T'y± SSC m c c 2 ~ 

15.6GeV L^ 4 6 ^; 1/4 (7^ SSC /350). After interacting with 
the MeV photons again, tertiary pairs will be produced 
with a typical Lorentz factor of ~ 180. They then re- 
scattcr various photon fields again. The detailed investi- 
gation of this cascade process has to be modeled numer- 
ically, and is beyond the scope of this paper. Nonethe- 
less, a generic picture is that at least 1/3 of the total 
energy of the reprocessed photons is radiated in the sub- 
GeV range, making these GRBs interesting targets for 
GLAST observations . 

1 In the Klein-Nishina regime, for one electron, the energy loss 
in one scattering is ~ 7 c m e c 2 , the energy loss rate is Pic ~ 
7 e m e c 2 A(a;)(T T n 7 c, where n 7 ~ U 1 /(e^ / /T) is the number den- 
sity of 7— ray photons. The synchrotron radiation power satisfies 
Psyn = (4/3)ctt7o/3 2 Ubc. The Compoton parameter can be esti- 
mated by Y KN = Pjc/Psyn ~ A(x)U-y / (xUb) , where we have taken 
7c = 7^! SSC . Since U^/Us is just the Compton parameter Y in the 
Thomson regime, we have Y KN ~ A(x)Y/x for x ^s> 1, i.e., in the 
extreme Klein-Nishina regime the IC cooling is very inefficient. 
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For typical parameters taken in this work (z = 1 and 
Q x = 1), the fluence of the sub-GeV flash from the very 
early FS can be estimated as 



2 

■2 x 1CT 6 ergs cm~ 2 , 



5~^e c r(r- l)m p c 2 (l + z) 



R x (wind) 



R 2 ndR/Dl 

(14) 

where is the luminosity distance, which as 2.2 x 
10 28 cm for z = 1. The fluence for the RS compo- 
nent may be comparable or less than this value. For 
EGRET, the estimated fluence threshold in the low in- 
tegration time (tint) regime (i;„t < 1.7 x 10 3 s) is ~ 2.1 x 
10~ 6 ergs cm~ 2 for a typical photon energy ~ 400 MeV 
(e.g. Zhang & Meszaros 2001). Our predicted fluence 
S ~ 10~ 6 ergs cm -2 at sub-GeV energy range is below 
the EGRET detection threshold, which explains the none 
detection of such flashes in the EGRET era. The fluence 
threshold of GLAST is roughly ~ 4x 10~ 7 ergs cm -2 for 
t int < 10 5 s (e.g. Zhang & Meszaros 2001). We there- 
fore expect that many sub-GeV photon flashes can be 
detected in the GLAST era. 

4. NEUTRINO EMISSION 

4.1. Neutrinos from photomeson interaction with the 
synchrotron photons in RS and FS 

The protons accelerated from both the FS and the RS 
would produce high energy neutrinos via photomeson in- 
teraction, mainly through A-resonance. In the RS re- 
gion, DL01 have found that the typical neutrino energy 
is 3 x 10 15 - 3 x 10 17 eV. 

We now estimate the neutrino emission from the FS 
region. As mentioned in ( see the paragraph below 
equation |1U|). the fresh electrons accelerated by the FS 
lose energy mainly by SSC cooling rather than by IC 
cooling. As a result, the flux above should be low- 
ered by a factor 1/(1 + Yssc)- The typical synchrotron 
radiation frequency (i^) of the FS reads 



3.1 x 10 18 Hz e 2 



£8.-2-^52.6^4* X R 1 $, (15) 



e,-l c B.-2 



The protons that interact with the hv m photons at the 
A— resonance have an observed energy of 



6^ obs (l)«0.3GeV 2 r 2 /[(l + z)^ ; 

rsI.8 x 10 17 cV (l + z)- 1 



e 2 e~ 1/2 

t e,-l t B,-2 



r -1/2 .1/2 „ 
-^52.6 A * it l 



(16) 

The photomeson interactions include (1) production of tt 
mesons, />7 — > p + ir° and jry — > n + 7r + , and (2) decay of 
tt mesons, ir° — > 27 and ir + — > + — > e + + v c + + 
Vfj,. When protons interact with the FS photons with 
typical energy hv^, these processes produce neutrinos 
with energy e^ obs (l) ~ O.Q5e^ obs (l) (e.g. Waxman & 
Bahcall 1997), i.e. 



< obs (l)~9x 10 15 eV (1 + z) 



-1-1/2 2 

t B,-2 t e,-: 



t -1/2 ,1/2 p 
^52.6 A * rt 15- 



(17) 

In principle, 7— rays of similar energies are produced by 
7T° decay. But as shown in the MeV photon flow 
is optically thick for these high energy photons, and a 
e^- pair cascade develops. The bulk 7— rays emerge at 



much lower energies ~ GeV, which adds to the sub-GeV 
emission discussed in 

Similar to Halzen & Hooper (2002) and Guetta et 
al.(2004), at e^. b s (l)j the fraction of energy loss of pro- 
tons through photomeson interaction with the hv l m pho- 
tons, /ti-(I), can be paramterized as (in the calculation, 



we have taken Yssc ~ 3.2e 



1/2 -1/2 

c,-l £ B,-2 



/ ff (l)~0'A n 7 (l)Ai?' < 

e_eE 

" 2(1 + Yssc) (hvtJAarRR; 

~n aa -3/2 r -1/2 .3/2 



-fTA < Xy_ 



> 



(18) 



where <ta ~ 5 x 10 _28 cm 2 is the cross section of the 
A— resonance, < x p ^ n >~ 0.2 is the average fraction 
of the energy transfered from the initial proton to the 
produced pion, AR' is the width of the shock, and 
n 7 (l) « (i?/i? x )[e ^/(2^ nl )]/[4^(l + Y SS c)i? 2 Ai?'] (for 
R < -Rx) is the photon number density, E ~ 4 x 10 53 ergs 
is the isotropic energy of the ejecta at the beginning of 
the afterglow, as defined in [J2 It is interesting to see 
that /tt(1) is independent on t since cx Br 1 , 

Below we estimate the maximum proton energy 
(e^ 1 ) accelerated by the shocks. In general, ep 4 = 
min[e^(l),e^(2),e^(3)] satisfies three constraints. (1) 
The comoving shock acceleration time t' a ~ e p /TeBc 
should be smaller than the comoving wind expansion 
time t^ ~ R/Tc, which yields e^(l) ~ eBR. The nu- 
merical value reads 



f (1) 



5 x 10 20 eV e^ 2 LliU* 1/4 - 



(19) 



In this work, the superscript ' represents the parameter 
measured in the comoving frame. 

(2) The comoving proton synchrotron cooling timescale 
^cooi = [67rmpC 3 /(TTTO 2 ]rep 1 B~ 2 should be longer than 
the comoving acceleration timescale t' a (e.g., Li et al. 
2002), which results in 



C(2) 



4 x 10 20 eV e-^Ll^A^Rli 2 . (20) 



(3) The comoving proton cooling timescale due to pho- 
tomeson interaction should also be longer than the 
comoving acceleration timescale t' a . t' w can be derived 
from f w ~ t'Jt'n (e.g. Waxman & Bachall 1997), i.e., 
f w ~ {R/Tc)f-\ where U = A(l)(e P /^) 1/2 for e p > 
el(l). Requiring t'„ = t' a yields e p /{TeBc) = R/[UTc], 
which can be simplified as £^(3) ~ f~ 1 e p A (l) ~ 
U(l)[e™(3)/e\;}-V 2 ef(l). We then have 

ef(3)^Miy 2/3 [ef(l)} 2/3 (^) 1/3 

c4xHPeV el^eZAiU^Rli 3 , (21) 

Since e£ )0bs < = e^(3) ~ 4 x 10 19 eV, the predicted 
flux of neutrino emission is insensitive to the actual value 

After the pions (muons) are generated, the high energy 
pions (muons) may lose energy via synchrotron emission 
before decaying, thus reducing the energy of the decay 
neutrinos (e.g. Guetta et al. 2004). For pions, this effect 
should be taken into account when the pion lifetime t' v ~ 
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loss time t' n syn 



6nm^c i /(aTm 2 e' 7T B 2 ). We can then 

/ 

7r,syn> 



define a critical pion energy e^. c by requiring = tl 
which gives 



2.6 x 10~ e'^/(m w c )s is comparable to the synchrotron where P^^^ ~ 6x 10~ (e^obs/3 x 10 eV) is the prob- 
ability that a neutrino produces a detectable high en- 
ergy muon for e„ i0b s > 10 3 TeV (Gaisser et al. 1995). 
The factor 1/8 takes into account the fact that there are 
three neutrinos produced with each carrying about 1/8 
of the energy loss of the proton, and that due to neutrino 
flavor oscillation, only 1/3 of the neutrinos are the de- 
tectable muon neutrinos. For a nearby GRB with z ~ 0.1 
(D L — 1.4 x 10 27 cm) and £~ 4 x 10 53 ergs, the result- 
ing detectable events is ~ 0.4. For such a nearby ener- 
getic GRB, the prompt 7— ray energy fluence would be 



1.1 x 10 17 eV e'^lL-^Al^R, 



(22) 



C B, -2^52. 6 

above which the neutrino emission is suppressed signifi- 
cantly. Correspondingly, we have a critical observed en- 
ergy for (since e' ^ w 4e'^), 

Te' c 



u u ,obs 



4(1 + z) 



1.2 x 10 i8 eV (- 



.1 



-1/2 



^15 



( 23 ) 

above which the slope of the spectrum steepens by 2 
since 4 iSyn /r4 oc e' n ~ 2 . 

Muons have a lifetime ~ 100 times longer than that of 
pions, so that the energy cutoff of and v c will therefore 
be 10 times smaller than e£ obs (e.g. Guetta et al. 2004). 

With equation (|23|) we have 



.Vf, I ,obs 



£ l/ u ,obs „ ml7 ir /^ Z \-l -l/2n 

1.2 x lO^eV (— ^— ) x e B L 2 R l5 , 



10 



(24) 

Equation Ijl8|l is valid only for protons at the break en- 
ergy £p obs (l). Considering emission from all the protons, 
the resulting neutrino spectrum would trace the photon 
spectrum. Generally one has 



/*(!)[; 



ii.i , , - b / 

(I) J ' fc p,obs ^ fc p ,obsV 

/-(l)[|^)] - 5 .<ob s (l) < £ P,obs < e° )0b (25) 

1 f° < f u < f M 



p,obs ^ C P>' 

/J> ^0.5 



p,obs' 



where e p,obs is defined as (e° >obs /e° i0bs ) 



VA(i); 



M _ M 
p,obs p 

2 d 2 N, 



Cp /(l + z). The neutrino spectrum finally reads 



'.obs 



1 (1 + z)L 
8 AnDl 



/7r(l)[ £ b b ,(l)] 1 ' 1 ' e iAobs < e ^obs(l)' 

/ir(l)[ e b ^ (!) ]°' 5 J e ^,obs(-'-) < e v,ohs < e ° i0 bs> 



1/, obs 



< £l/,obs < e° 

< £iAobs < e 



(26) 



ZAobs " 



where the factor 1/8 takes into account the fact that 
charged and neutral pions are produced with roughly 
equal probabilities, and that each neutrino carries ~ 1/4 

0.056° obs . In equation 

> M 



of the pion energy, i.e. e° obs 



5 v. obs 



'.obs 



the neutrino flux is for e„, bs > e !fobs- 

As in many publications, we assume that the protons in 
the FS distribute as dnjd^ v oc 7~ 2 for 7 p m < 7 p < 7 p ,m, 
where j p<m = (1 - e e )(r - l)/ln(7 PjM /7p, m ) and 7 PiM = 
e p M /rm p c 2 . For a neutrino detector with a surface area 
of Sdet ~ km 2 , the detectable events for one burst can 
be estimated by (taking R = 0.5i? x (wind)) 

AT m ~ 1 (1 ~ e ° )£; 1 S, P 
ln[4obs/<obs(l)]f-(l) + ln(<obs/<obs) 



ln(7 p , M /7p,m) 



0.002i? 5 3. 6j D L 2 8 34 , 



(27) 



JF 7 = (1 + z)E 1 /(AttD 2 ) ~ 10- 3 ergs cm^E^^D^ 15 . 
This means that only for those very bright GRBs can 
the predicted neutrinos be detected. Similar conclusions 
have been also drawn in Dermer & Atoyan (2003) and 
Guetta et al. (2004). We note that low redshift long 
GRBs are rare. So far we have detected at least three of 
them, but none of them is energetic enough for the above 
purpose. 

Another method to estimate the detectability of the 
neutrino early afterglow emission is to assume that GRBs 
are the dominant sources of high energy cosmic rays 
(from the internal or external shocks, Waxman 1995; Vi- 
etri 1995). This leads to a detectability about an order 
of magnitude higher (see e.g. DL01) than the value pre- 
sented above. One possible reason for the discrepancy 
is that this second method is based on the assumption 
that ultrahigh energy cosmic rays are mostly produced 
by GRBs, which might not be the case. 

The atmospheric neutrino background flux is (f>u,bkg ~ 
10 _12 (e l ,,obs/10 14 eV)" 5 / 2 cm- 2 s" 1 sr" 1 . The typical an- 
gular resolution of the planned neutrino telescopes at 
TeV energies is about 1 degree. For a typical GRB 
with T 90 ~ 20s, £ ~ 4 x 10 53 ergs and z = 1, at 
e<y,obs ~5x 10 15 eV, the background neutrino number in 
this time span is ~ 5 x 10" 9 , which is < 0.002 (eq.[57|). 
If a ~ 5 x 10 15 eV neutrino well correlated with one GRB 
both in time and coordinates has been detected, the 
chance probability of being an atomospheric neutrino is 
only ~ 3 x 10~ 5 . Therefore, these GRB neutrinos can 
be detected above the background (see also Meszaros & 
Waxman 2001). 

4.2. Neutrinos from photomeson interaction with the 
MeV photon flow: ~ lO^eV 

Because of the significant overlapping of the MeV pho- 
ton flow and the shocked regions, high energy protons 
with energy larger than 2 x 10 15 eV accelerated from both 
the RS and the FS would interact with the MeV flow ef- 
fectively. We now study the neutrino emission process in 
this case. 

In the rest frame of the shocked regions, 7— ray photons 
with the observed typical energy obs w 300keV have 

an energy ejj rts (l + z)e^ obs /r. The protons that interact 
with the photons with energy at the A-resonance have 
an observed energy of 



eb obs (2)«0.3GeV 2 r 2 /(l + z) 2 e !; iObs 



2 x 10 15 eV Ll f 2 2 6 A* 



l/2^ 1+z -p2^ e y,abe 



300koV 



r\ (28) 



This energy is well below the observed maximum proton 
energy ep 4 (equations E]-jU). The generated neutrinos 
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have a typical energy 



^i/jObs 



(2) 



10 14 eV Ll^A 



"7, obs \_ 



\29) 



2 ' v 300kcV y 

Similar to H4.ll at e£ b s (2), the fraction of the energy 
loss of protons through photomeson interaction with the 
MeV photons can be parameterized as 



U{2) w 0.18L 7 , 52j L 



1/2.1/2 



52.6 



(1 + *)<o b B 

600keV 



The MeV flow photons has a broken power law spectrum, 
i.e., n(e 7 ) oc e" 1 below e 7! ob s and n(e 7 ) oc e~ 2 above 
e 7 obs- ^° ^ ne g enera l form of the fraction of the energy 
loss of protons through photomeson interaction with the 
MeV photons is 



h « A (2) 



1. 



a (2) 



£p,obs > e p^obs(2)j 
£p,obs < e p,obs(2) 



(31) 



As a result, the differential neutrino spectrum is also 
a broken power law, i.e. n(e v ) oc e^ 1 below e^ obs (2), 

and rt(e„) oc e~ 2 above e^ obs (2) (see also Waxman & 
Bachall 1997). However, for e^obs > e ° bs> the energy 
of the corresponding protons has been lost mainly by 
interacting with the keV soft photons from the FS region 
(see M.ll for detail). So the neutrino spectrum can not 
extend to e„ i0b s > e° obs - 

Similar to equation l|27|l . we now estimate the de- 
tectability of such 10 14 eV neutrinos in the early afterglow 
phase. For one typical GRB, taking R = 0.5i? x (wind), 
one has 



V CV cnts(2) ~ ~~ 



8 <obs(2) ^Dt 



2 ^dctPy — ■ fi 



ln[< obs /<obs(2)]/,(2) 

hl(7p,M/7p,m) 



0.001£ B3 . B £> 



-2 



53.6-L-Y 28.34 



, (32) 



where P v ^ ~ 1.7 x lO^e^obs/lTeV) 08 (e.g., Halzen & 
Hooper 2002) for ITeV < e„' obs < 10 3 TeV. The resulting 
detectability is comparable with that of 5 x 10 15 eV (see 
equation |27j). These neutrinos can be only marginally 
detectable above the atmospheric neutrino background. 
For a typical GRB with T 90 ~ 20s, E~ 4 x 10 53 ergs and 
z = 1, the 10 14 eV neutrino background event number 
is - 9 x 10~ 5 in the time span. It is < 0.001 (ea.|32]). 
but not significantly. For a detected ~ 10 14 eV neutrino 
well correlated with one GRB, the chance probability of 
being an atomosphcric neutrino is ~ 10%. 

The total detectable events for one typical GRB born 
in the wind environment (z — 1) is ~ 0.003, this rate 
can be also enhanced significantly when fluctuations in 
distance, energy are considered (e.g., Haltzen & Hooper 
2002). Therefore, nearby GRBs such as GRB 030329 are 
of great interest (e.g., Razzaque et al. 2004). Thanks 
to the correlation of the neutrinos with GRBs, both in 
time and coordinates, these GRB neutrinos are likely to 
be detected with the planned km 3 high energy neutrino 
detectors such as ICECUBE, ANTARES, and NESTOR. 

One caveat is that the ~ 10 14 eV neutrinos would also 
overlap with the ~5x 10 14 eV neutrinos generated from 
the internal shocks (e.g., Waxman & Bachall 1997), so 



that it might not be easy to distinguish both species. We 
notice, however, that if GRB ejecta are strongly mag- 
netized so that the prompt emission is not powered by 
internal shocks, the 10 14 eV neutrinos discussed here are 
still produced. 

5. DISCUSSION AND CONCLUSION 

Massive stars are likely the progenitor of long, soft 
GRBs. For an ultra-relativistic fireball emerging from 
a massive star and expanding into a dense stellar wind 
medium, a strong relativistic RS develops so early that 
both the RS and FS regions overlap in time and space 
with the prompt MeV photon flow. This overlapping 
leads to modification of the early RS and FS emission 
significantly, since the dominant cooling process for elec- 
trons is likely the IC with the MeV photons (especially 
for the RS region, Bcloborodov 2005; see also Fan et al. 
2005a for a more detailed calculation on the modification 
of very early R-band lightcurve) . 

Due to the tight overlapping of the MeV photon flow 
and the shocked regions (both FS and RS), the opti- 
cal depth of the multi-GeV photons is larger than unity. 
These high energy photons (SSC and/or IC emission 
components from the FS and RS regions) are absorbed by 
the MeV photon flow and therefore generate relativistic 
pairs. These pairs re-scatter the soft X-rays in the 
FS region and the prompt 7— rays and power detectable 
high energy emission. A significant portion of the re- 
processed photon energy is distributed in the sub-GeV 
band, making wind-interaction GRBs interesting targets 
for the GLAST. In the wind model, the total thermal 
energy contained in the FS region is larger or at least 
comparable with that contained in RS region. The above 
effect is therefore generic, regardless of the possible mag- 
netization of the ejcta (e.g. Fan et al 2002; Zhang et al 
2003, in which case the IC emission as well as the SSC 
emission of RS are suppressed). The predicted fluence of 
these sub-GeV flashes is high enough to be detected by 
GLAST. 

Equally interestingly, the early photon-shock interac- 
tion leads to interesting neutrino emission signatures. 
Besides the conventional neutrino component powered 
by the photomeson interaction in the RS region (DL01), 
in this paper we studied two new neutrino emission com- 
ponents, i.e. a~5x 10 15 eV component powered by the 
photomeson interaction in the FS region, and a ~ 10 14 
eV component powered by photomeson interaction with 
the MeV photon flow. The detectability of the two com- 
ponents is comparable. The spectrum of 10 14 eV neutrino 
is similar to that of the neutrinos generated in the inter- 
nal shock phase since the seed photons are nearly the 
same, although the typical observed energy is softened 
by a factor ~ (r/r ) 2 . The spectrum of the 5 x 10 15 eV 
neutrinos takes a more complicated form (see equation 
1261*1 . A caveat is that these neutrinos overlap with the 
~ 5 x 10 eV internal shock neutrino-component (Wax- 
man & Bahcall 1997), and hence, is difficult to distin- 
guish. On the other hand, they are generic regardless 
of whether the MeV emission is powered by the internal 
shocks or not. 

In this work, it is assumed that the wind medium is 
at rest at any radius. This may not be the case if the 
radiation front of the GRB is taken into account. As 
shown in Beloborodov (2002), due to the large amount of 



High energy photon and neutrino emission from GRBs 



7 



e pair loading (The pairs are created by the interaction 
of the prompt 7— rays with the back-scattered 7-rays by 
the wind medium), the medium will be accelerated and 
the ejecta moves in a cavity until it reaches a radius 
R gap ~ 3.3 x 10 15 £^ 5 2 3 (e.g., Beloborodov 2002). As a 
result, the shock crossing radius (see equation J3Jl) should 
be ~ -Rgap + i? x (wind) ~ 2i? x (wind). The overlapping 
of the MeV photon flow and the shocked regions is still 
significant since the ejecta essentially does not decelerate 
at R < Rgap. However, f v (l) and f v (2) involoved in the 
above calculation should be smaller by a factor of about 
1/4, so are the predicted neutrino fluxes. The optical 
depth for the multi-GeV photons (see eo.|12|) is ~ 1, so 
that they are still absorbed by the MeV photon flow and 
generate relativistic pairs, which re-scatter the initial 
MeV 7— rays to the sub-GeV energy range. Notice that 
the introduction of the radiation front will soften the FS 
shock emission significantly since the upstream of the FS 
is rich. If we assume that the number of the pairs is 
k times that of the electrons associated with protons, 7, f n 
(see eq.|S]) should be reduced to 7^/(1 + 2k). If k is as 
large as several tens, the elecrons accelerated in the FS 
will be cooled mainly via IC with the MeV photon flow 
rather than synchrotron radiation and SSC emission. A 
significant part of that IC component may be in sub- 
GeV energy range. Therefore sub-GeV flashes are still 
expected when the radiation front effect is taken into 
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account. 

In this work, only the wind model is discussed. In the 
ISM model, if the GRB is long enough, the overlapping 
of the MeV photon flow and the shocked regions could be 
also tight. In this case, R x (ISM) ~ 10 17 cm, and the opti- 
cal depth of GeV photons is very small (since r 77 oc i?~ 2 , 
see eq.^2])- I n such a case, GeV-TeV flashes are ex- 
pected, as suggested by Beloborodov (2005). Therefore, 
the sub-GeV photon flashes are a signature for GRBs 
born in a stellar wind. 

We note that the overlapping effect increases the neu- 
trino flux predicted from the collapsar GRB model to 
be close to the one predicted in the supranova model (cf. 
Dermer & Atoyan 2003). This is because the overlapping 
MeV photons effectively provide a convenient target for 
photomeson interactions, similar to the external photons 
involved in the supranova model. 
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